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BETTER BOUNDS FOR PLANAR SETS AVOIDING UNIT 

DISTANCES 

TAMAS KELETI, MATE MATOLCSI, FERNANDO MARIO DE OLIVEIRA FILHO, 

AND IMRE Z. RUZSA 


Abstract. A 1-avoiding set is a subset of W that does not contain pairs of 
points at distance 1. Let denote the maximum fraction of R n that 

can be covered by a measurable 1-avoiding set. We prove two results. First, 
we show that any 1-avoiding set in R n (n > 2) that displays block structure 
(i.e., is made up of blocks such that the distance between any two points from 
the same block is less than 1 and points from distinct blocks lie farther than 1 
unit of distance apart from each other) has density strictly less than l/2 n . 
For the special case of sets with block structure this proves a conjecture of 
Erdos asserting that mi(R 2 ) <1/4. Second, we use linear programming and 
harmonic analysis to show that mi(l 2 ) < 0.258795. 


1. Introduction 


The unit-distance graph of R” is the graph whose vertex set is R ra and in which x 
and y are adjacent if ||cc — j/|| =1. A well-known problem in geometry, going back to 
Nelson and Hadwiger (see Soifer [T2] for a historical survey), asks for the chromatic 
number %(R") of the unit-distance graph of R n . 

A related problem considers independent sets of the unit-distance graph. Let G = 
(V, E) be a graph. A set I C V is independent if it does not contain a pair of ad¬ 
jacent vertices. The independence number of G, denoted by a(G), is the maximal 
cardinality of an independent set. 

A set A C R" is independent in the unit-distance graph if it does not contain 
pairs of points at distance 1, that is, ||cc — y\\ ^ 1 for all x, y € A. We also say that A 
avoids distance 1 or that it is a 1-avoiding set. A possible measure for the size of 
an independent set in this case is its density, that is, the fraction of space that it 


covers (see 11.1 for a rigorous definition). Our aim is to estimate the maximal (or 
more precisely the least upper bound of the) fraction of R n that can be covered by 
a measurable set that avoids distance 1. 

We denote this maximal fraction by mi(R n ). In ^3j we show that mi(R 2 ) < 
0.258795. This result is related to a conjecture of Erdos (cf. Szekely |Hj), stating 
that toi(R 2 ) < 1/4. 

As for the chromatic number of the Euclidean plane, all that is known is that 4 < 
y(R 2 ) < 7. The upper bound comes from a simple periodic coloring of R 2 , whereas 
the lower bound comes from a finite subgraph of the unit-distance graph, the Moser 
spindle (cf. Moser and Moser 0), whose chromatic number is 4 (see Figure [l]). 
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Figure 1. On the left, the Moser spindle j5j. Each segment has 
length exactly 1. On the right, the optimal tortoise in Croft’s 
construction [2]. 


In view of the difficulty of computing \(R ra ), Falconer [3] introduced the measur¬ 
able chromatic number , denoted by x m (R"), in which the restriction is added that 
the color classes must be Lebesgue-measurable sets. In other words, one wishes to 
partition M" into the minimum possible number of Lebesgue-measurable 1-avoiding 
sets. Obviously, x m (M n ) > x(R"). Falconer proved that x m (M 1 2 ) > 5. Since 

m 1 (R") Xm (M") > 1, 

showing Erdos’ conjecture would give another proof of Falconer’s result. 

A simple lower bound for mi(R 2 ) comes from the following construction. Con¬ 
sider the hexagonal lattice with minimal vectors of length 2 and place at each 
point of the lattice an open disk of radius 1/2. This is a 1-avoiding set of den¬ 
sity 7r/(8 v/ 3) = 0.2267.... A slight improvement was given by Croft [2]. His 
construction is as follows. Shrink the hexagonal lattice slightly so as to have min¬ 
imal vectors of length 1 + x, where x < 1, and place at each lattice point the 
intersection of an open disk of radius 1/2 and an open regular hexagon of height x. 
The disks guarantee that inside each block every distance is less than 1, while the 
hexagons guarantee that points from different blocks have distance greater than 1. 
Taking x = 0.96553 ... maximizes the density of the union of all the blocks, showing 
that mi(M 2 ) > 0.22936. This intersection of a disk with a hexagon has often been 
called a tortoise ; Figure |T| shows an optimal tortoise. 

Any finite subgraph G = (V,E) of the unit-distance graph of M n provides an 
upper bound for mi(R n ), namely a(G)/ |Vj; this can be seen via a simple averaging 
argument^ For the plane, we then immediately have that mi(R 2 ) < 1/3, as can 
be seen from the equilateral triangle. A better bound of 2/7 is provided by the 
Moser spindle; this is the best upper bound that has been obtained from a finite 
subgraph. 

Using further ideas, Szekely 13] proved a bound of « 0.279 < 2/7 « 0.285. 
Oliveira and Vallentin HD] gave the currently best known upper bound of ~ 0.268. 
Their method is based on a mix of linear programming and harmonic analysis; it 
is a strengthening of it that will be used in fj3] to prove that toi(R 2 ) < 0.258795. 

Frankl and Wilson J4] construct finite subgraphs of the unit-distance graph of R n 
whose chromatic numbers grow exponentially fast in n. These same subgraphs can 
be used to provide the asymptotic upper bound 

mi(R n ) < (1 + o(l))1.207 _ra . 


1 I bis is related to the following observation: Let G be a subgraph of a finite vertex-transitive 

graph H. Then a{H)/\V(H)\ < a{G)/\V(G)\. 










Better bounds for planar sets avoiding unit distances 


3 


via the averaging argument. The method of Oliveira and Vallentin uni also provides 
an exponential upper bound, mi(R") < (1 + o(l))1.1654 _ ", which is somewhat 
weaker than the bound of Frankl and Wilson. A combination of both arguments, 
requiring detailed analysis, was used by Bachoc, Passuello, and Thiery [I] to obtain 
the best known asymptotic upper bound 


mi(K”) < (1 + o(l))1.268“". 

We remark here that the results of [Tj and the present paper do not overlap: the 
authors of QQ obtain asymptotic upper bounds as n —> oo, and some numerical 
improvements for dimensions 4 < n < 24. 

The behavior of x(R") under restrictions placed on the color classes has also 
been studied. We have already mentioned Falconer’s measurable chromatic number. 
Other restrictions include requiring all classes to be either open or closed sets — 
in both cases it can be shown that the chromatic number of R 2 is either 6 or 7 
(cf. Soifer [12]). 

Similarly, we may place some natural restrictions on 1-avoiding sets and study 
their maximal possible upper densities. In particular, we will consider sets with 
block structure. We say that a set A has block structure if it is a union 

OO 

A= {JA Z 

i—0 


of blocks Ai, where ||x — y\\ < 1 if x and y belong to the same block, and ||x —y|| > 1 
if x and y belong to different blocks. Notice that we do not require A to be 
measurable (but this will not make an essential difference; cf. the proof of Theorem 


2.1 below). 


A set with block structure is clearly a 1-avoiding set. All known constructions 
of 1-avoiding sets of “high density” are actually constructions of sets with block 
structure, like the hexagonal lattice construction of disks, or Croft’s construction. 
Recall that Erdos conjectured that toi(R 2 ) <1/4. Larman and Rogers, and before 
them Moser (cf. Larman and Rogers [7]), made the following conjecture: the volume 
of a closed 1-avoiding set inside a ball of radius 1 in R" is less than 1/2" of the 
volume of the ball. A simple argument shows that this conjecture implies mi(R") < 
1/2", and it is therefore a generalization of Erdos’ conjecture. In (J2]we will show 
that any subset of R" (n > 2) with block structure has upper density less than 1/2". 


1.1. Preliminaries and notation. Throughout the paper, A(A) and A(.4) will 
denote the Lebesgue measure of a (measurable) set A, and the outer measure of 
any set A, respectively. 

Let A C R" be a measurable set. We say that its density is 6(A) if for all p £ R" 
we have 


<5(A) = 


lim 

r—> oo 


A(An S(p,r)) 


\(S(p,r)) 

where S(p 7 r) is the n-dimensional cube of side 2r centered at p. For a set that has 
a density, the cube can be substituted by any reasonable body, like the ball, say, 
without changing the resulting density. 

A measurable set may not have a density, but every set has an upper density 


We define 


5(A) = lim sup 

r—> oo 


X(A n S(p,r)) 
A (S(p,r)) 


mi(R") = sup{ 5(A) : A C R" is 1-avoiding and measurable}. 

We say that a set A is periodic if there is a lattice L C R n that leaves A invariant, 
that is, x + A = A for all x £ L. Then L is the periodicity lattice of A. 
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Measurable periodic sets have densities. Moreover, a simple argument shows 
that the densities of periodic 1-avoiding sets can come as close as desired to mi(M”) 
(cf. Oliveira and Vallentin |10jb So when computing upper bounds for mi(R”) we 
may restrict ourselves to periodic sets. 

2. Sets with block structure 

In any dimension, all known examples of 1-avoiding sets of “high density” are 
made up of disjoint blocks, i.e., they are sets of block structure as defined in the 
introduction (see the end of this section for an example of a 1-avoiding set of small 
positive density which does not have block structure). 

On the one hand it is natural to consider this class of sets because human imag¬ 
ination of 1-avoiding sets seems to be more or less restricted to this class. On 
the other hand, it seems very elusive to prove rigorously that a 1-avoiding set of 
maximum or close to maximum density must be block-structured. 

The following theorem, for n = 2, is a special case of the conjecture of Erdos 
(cf. Szekely [2]) given in the introduction; for n > 3, it is a special case of a con¬ 
jecture of Larman, Rogers, and Moser (cf. Larman and Rogers [7]), also discussed 
in the introduction. 

Theorem 2.1. Let n > 2 and let A C R" be a 1-avoiding set having block structure. 
Then 5(A) < 1/2” — e n for some e n > 0. 

We will prove the following slightly stronger result. 

Theorem 2.2. Let n > 2 and let A\, A 2 , ■.. C R" be sets of diameter at most 1 
such that the distance of any two of them is at least 1. Then the upper density of 
A = (J“i Ai is at most 1/2" — e n for some e n > 0. 

Observe that both properties (diameter at most 1, distance at least 1) are pre¬ 
served if we replace the sets by their closure, so we may assume that the sets Ai 
are all closed. This also makes them measurable. (This reduction works only under 
the assumption of block structure. In the general case, it is fairly easy to show the 
existence of a (non-measurable) 1-avoiding set of full outer measure, that is, such 
that the inner measure of its complement is 0.) 

First we present a simple proof for the weaker statement when e n is discarded, 
that is, 5(A) < 1/2”. 

Let B r denote the open ball of radius r around the origin and let 
Ci = Ai + U 1/2 = {a + b : a £ Ai, b £ B 1/2 }. 

It is clear from the assumptions that CiC\Cj =0, for all i ^ j. 

Applying the Brunn-Minkowski inequality (see e.g. [51 Theorem 4.1]) to the sets 
Ai, B 1 /2 we obtain A(Cj) 1 /" > A(A i ) 1 /”+A(I? 1 / 2 ) 1 ^ n - Furthermore the isodiametric 
inequality gives A (Ai) < X(B 1 / 2 ). By combining these inequalities we get 

A(Aj) 1 /” A (A,) 1 /” = 1 1 _ 1 

J A(C i ) 1 /” - A(Aj) 1 /” + A(B 1/2 ) 1 /” (\{B 1/2 )\ x/n 1 + 1 2' 

1 + ^ A (A t ) ) 

Since the sets Ci are pairwise disjoint, Ai C Ci, and A = (J !; Aj, this shows that 
5(A) < 1/2”. 

The plan to show 5(A) < 1/2" — e„ is the following: If 5(A) is close to 1/2” then 
in the above argument we must have that for most i the isodiametric inequality is 
almost an equality and diamA^ is close to 1. By a stability theorem this implies 
that each such Aj is very close to a ball of radius 1/2 and then most of the sets 
Ci = Ai+Bi /2 are very close to unit balls. But the density of any unit ball packing 
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is well-separated from 1, and then so is the density of (J ?; C*. Since A has density at 
most 1/2" in Cj, this implies that the density of A is well separated from 1/2". 

The stability result we use is the following theorem of Maggi, Ponsiglione, and 
Pratelli [3]. 


Lemma 2.3. Let E C R" be a measurable set with X(E) > 0 and diami? = 2. 
Then there exist x, y £ R" such that 


E C B(x, 1 + r) and B(y, 1) C E + B r , 
where B(z,R) denotes the ball centered at z with radius R and 


r = K, 


(KBi 


V A (E) 

for some constant K n that depends only on n. 


- 1 


l/r, 


Note that for sets of diameter 2 the expression X(Bi)/X(E) — 1, which is called 
isodiametric deficit by Maggi, Ponsiglione, and Pratelli, is nonnegative by the iso- 
diametric inequality and expresses the error in that inequality. 

We will need the following simple corollary of the above result. 


Corollary 2.4. For any n> 2 there exists an increasing function f3 = /3 n : (0, oo) —► 
(0, (X)) with limp^o /3 n (p) = 0 with the following property. For every measurable 
E C R" with X (E) > 0 and diam I? < 1 there exist x, y £ R" such that 


E C B(x, 1/2 +/3(p(E))) and B{y,1/2) C E + B p{p{E)) , 
P{E) = ~ L 


where 


Proof. By rescaling Lemma |2.3| we get that for any measurable E C R" with 
A (E) > 0 and diami? < oo there exist x, y £ R" such that 


(2) ECBlx, 


where 


diami? 


+ r I and B ( y, 


diam E \ 


2 ) 


E. E -\- B r , 


= K n 


diami? (X (B( diamS )/ 2 ) 

~Ee) 


- 1 


1/n 


If diami? < 1 then 


(3) 

SO 

(4) 


= ( A(B <" ) -' 2) - l) 1 '" < 


E C B \ x, 


diam E 


1 K 


Fntjr: i 1 /" 


+ r) CB[x,- + —p(E) 


By the isodiametric inequality we have A (E) < X(B(dieimE)/ 2 )i so 

1 


o(E) = X<yBl/2 ' > - 1 > X ( Bl / 2 ^ 

X(E) ~ A(B( diamB) / 2 ) * (diam E) r 


- 1 = 


- 1 , 


hence 

<5) diamE - (ryh) 

From the second part of using B(z, a) + B b = B(z, a + b) we get that 
(6) B(y,l/2)CE + B rHl — diam E) / 2- 
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Then 0,0.0 and ([6]) show that the function 


Mp) = ^ + Ui- 


1 + p 


l/r 


has all the required properties. 
Now we prove Theorem 1 2.2 1 


□ 


Proof of Theorem \2.2\ We fix the dimension n; all constants may depend on it. We 
will show that if N is large enough then the density of A in [-N, N] n is at most 
1/2" — e n for some positive e n . 

Suppose without loss of generality that the blocks Ai are enumerated so that 
A-[. ..., A rn are the ones that have nonempty intersection with [-N, N] n and so 


(7) 

Let 


A n [-N, N] n = (J (A t n [-N, TV]"). 


i= 1 


, 4 K B l/2) 1 
Pi = p(Ai) = . , . . -1. 


Then using 0 we get 
( 8 ) A ^ V " 


A (Ai) 


A(C' i ) 1 A 


< 


1 


1 


1 + 


(x (Bi,A\ 1/n l + (p l + l) 1 /" 

V HAi) 


Let a and p be two small positive constants that will be specified later. Let 
I = {i £ {1,... ,n} : pi> p} and J = {1,..., n} \ I. 

By 0, for every jg / we have 

A(Tb) < A (Ci) ( 1 + (/9+ iji/n) • 

By 0. for every i £ J we have A(Aj) < A(Ci)/2 n . Therefore we have 

n [-"• "1”) < E m) ( 1 + ( / +1) v ) " + E^ 


ieJ 


sE^-E*«y(£-( T 


i =1 i£l 

First we consider the case when 


+ (p + 1)!A 


^A(Ci) > a{2N) n . 


iei 


Using that the sets Ci are pairwise disjoint and (J™ x Q £ [—N — 2, A/" + 2] n , we get 

(9) \{An[-N,N] n ) < (2(iV + 2)r 1 a 


2=1 

A([—TV, N] n ) ~ (2 N) r ni ' 

Next consider the case when 

(io) E A (Ci) < a(2N) n . 

ier 


l + O + l) 1 /!; 


By definition for any i £ J we have p(Ai) < p. By Corollary 2.4 this implies that 
for any i £ J there exist Xi and yi such that 

(11) Ai C B(xi, 1/2 + /3(p)) and B{yi, 1/2) C + B^. 

Since linip^o (3(p) =0 we can guarantee f3(p) < 1/2 by taking p small enough. 
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Note that the first part of (11) implies 

(12) G,=i,+5 1/2 C%1+^)) 
and the second part of 0 implies 

(13) B(y i ,l-0(p))QA i + B 1/2 =C i . 

Since the sets Ci are pairwise disjoint, this implies that the balls B{y il 1 — f3(p)), 
for i £ J, are pairwise disjoint. 

Let A n be an upper bound on the density of the union of disjoint balls of the 
same size in R n . Although the best A„ is not known for general n, it is known that 
A„ < 1 for n > 2. 

Thus the density of (J i£j B{y i ,l~ j3(p )) in [ —N — 2, N + 2] n is at most A n , that 
is, 

(14) ( 2(JV + 2))- ^ A(Blil ‘’ 1 ~ W)) - A "- 

Notice that X(B(xi,l + P(p)))/X{B(y u 1 - /3(p))) = ((1 + /3(p))/{ 1 - /3(p))) n . 
Let Aj = (J iG jAi. Then, using first <[T|) then 0 and finally fl4|) we get that 


A(A/) 


< 


A(C'i) < 


1 


X([-N,N] n ) - (2N) n ^ 2 
(15) ieJ 


(2 Ny 


E 


X(B(xi,l + fi{p))) 


< 


i£j 

(2(N + 2)) r 


(2 Ny 


•A, 


i +m y ± 

i -m) 2n ' 


Let Ai = U ieJ Ai. By (jTj) and (10) we have 
(16) 


X(Ai) 


a 

X{[-N,N] n ) - 2™’ 


Combining (|15|) and (|16|) we get in this case 
(17) 


X{An[-N,N] n ) ^ (2(7V + 2)) n A 

^ - -- ’ 


1 +^(p) V . , a_ 

l-P(p)J ’ 2 n + 2 n ' 


A {[-N : N] n ) ~ (2 NY 

Finally, choose the positive constants a,p and s n so that j3(p) < 1/2, 

i+MV 


A, 


and 


i -m) 


i 


< a I — - 


T ol 2 n £ n < 1 


1 


i + (p + iyA 


Then by (|9j) and (17) we get that <5(A) < l/2 n — e n in both cases, which completes 
the proof. □ 


The theorem above suggests that one should try to prove that any 1-avoiding 
set of “high density” must have block structure. A natural idea is to take any 
1-avoiding set A and try to modify it in some manner to obtain a new 1-avoiding 
set A having block structure and at least the same density as A. Unfortunately, we 
could not prove anything rigorous along these lines. 

We end this section by presenting an example of a 1-avoiding set of positive (but 
small) density which does not have block-structure. This example also shows that 
not every 1-avoiding set of positive density can be modified in a natural way to 
obtain a new 1-avoiding set that has block structure and larger or equal density. 

Consider the scaled integer lattice (cZ) 2 C M 2 with c = 2\/2 — 2, and place 
an open disk of radius r = (3 — 2y/2)/2 at each lattice point. It is easy to check 
that the distance between points of disks around adjacent lattice points is less than 
1 and the distance between points of disks around nonadjacent lattice points is 
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bigger than 1. Therefore this is a 1-avoiding set without block structure. Simple 
calculation shows that the the density of this example is <5 = r 2 7r/c 2 « 0.0337. (A 
somewhat higher density could be achieved with a similar construction using the 
hexagonal lattice instead of the integer lattice.) 

3. A BETTER UPPER BOUND IN THE PLANE 

We now show how a strengthening of the method of Oliveira and Vallentin [TT)1 
can be used to provide the bound below. 

Theorem 3.1. 

TOi(R 2 ) < 0.258795. 

Let A C R" be a measurable and periodic 1-avoiding set with periodicity lat¬ 
tice L C R n . Its autocorrelation function is the function /: R" —> M defined by 

f(x) = S(A n (A — x)). 

Determining mi(R n ) is equivalent to the following optimization problem: find 
a function / that maximizes /(0) and is the autocorrelation function of a periodic 
1-avoiding set. The difficulty here lies in the fact that we do not have a charac¬ 
terization of autocorrelation functions of 1-avoiding sets. If we give up on finding 
autocorrelation functions, but settle for functions satisfying a few of the constraints 
that autocorrelation functions do, then we get a relaxation of our original problem 
and an upper bound for mi(R n ). The following lemma gives some such constraints. 
Recall that by a(G) we denote the independence number of a graph G. 

Lemma 3.2. Let f be the autocorrelation function of a measurable and periodic 
1-avoiding set A C R". Then: 

(1) f(x) = 0 if\\x\\ = 1; 

(2) if G = (V, E) is a finite, nonempty subgraph of the unit-distance graph 
of R n , then 

E /(*) < /(O)or(G); 

(3) if C C R n is a finite set of points, then 

E f(*-v)>\c\m-i, 

f*>i/}e(?) 

where (^) is the set of all pairs of points in C. 

Proof. For (1) it suffices to observe that, since A is 1-avoiding, if ||x|| = 1 then 

A n (A — x) = 0. 

For (2), we claim that any 2 £ R" belongs to at most a(G) of the sets A — x 
for x £ V. Indeed, say z belongs to all sets A — x t for {x\, ... ,Xk} Q V. This 
means that cu — Xi = z with some cu £ A. Now if k > a(G'), then there is a pair of 
points Xi, Xj adjacent in G , which means \\ai — a^H = ||a:* — cc_j || = 1, a contradiction. 
This observation now gives 

6(A) > (J (A n (A — x)) J > a(G)- 1 E (5 (^ n ( 2l - ;c )). 

EgF ' x£V 


as wanted. 
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Property (3) is an application of the inclusion-exclusion principle. We have 

1 > sf |J A-. 

\ec 


6((A - x) n (A - y)) 

xeG Ks/le(?) 

= \C\S(A)- £ S (An(A-(x-y))), 

U,y}e(£) 

and together with the definition of / we are done. 


□ 


Constraint (2) was observed by Oliveira and Vallentin pj] in the case when the 
graph G is a clique. Constraint (3) was used by Szekely [131 . 

Remark 3.3. The autocorrelation function encodes information about pairs of 
points in the set. We may also consider higher-order correlation functions which 
encode information about tuples of points. 

For every k > 0, consider the function F k : (R™) fe x {0, l} fc+1 —► R given by 

Fk{{x i,. ..Xk), (eo,ei,.. .,£*,)) = 5{A e ° n (A - x i) £1 D • • • D (A - x k ) ek ), 

where the notation X 1 = X and X° = R" \ X is used. 

These functions are nonnegative and satisfy the following self-consistency rela¬ 
tion: for all k > 1, xi, ..., x k and £q, E k , ..., £ k -i we have 

(18) F k —\ ((*^1 ) * * * 5 X k — 1 ) , (^Q , . . . , Ek— 1 ) ) 

= F k ((x 1 ,.. .,x k -i,Xk), (eo, ■ ■. ,£fc-i,0)) 

+ F k ((x i,.. .,x k -i,x k ), (eo, • • • ,£fc- 1 , !))■ 
Also, similarly to property (1) in the lemma above, we have 
Fk{(x i,.. -,x k ), (e 0 ,£i,.. .,£*,)) = 0 

whenever there exist indices i, j such that ||Xi — Xj || = 1 and Ei = Ej = 1. This 
property and the self-consistency relations together imply properties (2) and (3) 
in the lemma above (we omit the proof of this fact). Therefore, one may regard 
the self-consistency of higher-order correlation functions as the common source of 
properties (2) and (3). o 

Remark 3.4. Using Kai Lai Chung’s inequalities as Szekely and Wormald did in 
|16j . one can obtain a generalization of property (3) in Lemma |3.2| Namely, for 
every m > 1 we have j / } g ( c ) f( x — y) > m|C|/(0) — (and, in particular, 

to = 1 yields property (3)). In principle, these inequalities could be used to further 
improve the upper bound on TOi(R 2 ), as was the case in higher dimensions in [16] . 
However, when we implemented the inequalities in the computer code (cf. the 


description of section 3.1 below), they did not yield any numerical improvement, o 


In order to construct an optimization problem in which an autocorrelation func¬ 
tion of a 1-avoiding set is to be found, we shall parametrize such functions via their 
Fourier series. This parametrization will also suggest one further constraint satis¬ 
fied by autocorrelation functions. All the facts we state from harmonic analysis are 
quite standard; the reader looking for reference is advised to consult the book by 
Katznelson [B]. 

Given measurable functions /, g: R" -4 C, write 

(!9) (f,g)= J™ / f(x)g(x)dx , 


T—foo (2 T) n J[ 


— T,T1’ 
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when the limit exists. 

Let L C R n be a lattice. A function /: R™ —>• C is periodic, if it is invariant 
under the action of L: we have f(x + v) = /( x) for all x £ R” and v £ L. Lattice L 
is a periodicity lattice of /. Such functions are in a natural 1-1 correspondence with 
functions f: W 1 /L —> C; with an abuse of notation we will use the same letter for 
both. 


Now (19) defines an inner product in the space of square-integrable, complex¬ 
valued functions with periodicity lattice L, isomorphic to L 2 (R n /L). Equipped 
with this inner product, L 2 (R n /L) is a Hilbert space. Functions 


where L* = {u £ R ra : u ■ v £ Z for all v £ L } is the dual lattice of L, form a 
complete orthonormal system of L 2 (R n /L). 

The Fourier coefficient of / at u £ 2nL* is f(u) = ( f,\u )■ Since the \u form a 
complete orthonormal system, we have the Fourier inversion formula 

/(*)= E hu)e iu ' x , 

u£_2tvL* 


with L 2 convergence. 

Let A C R n be a measurable set with periodicity lattice L. We denote by 1 a 
the indicator function of A, an L-periodic function. The autocorrelation function 
of A is 


f(x) = ( 1 A , 1 A-x)- 


Since /( x) is expressed in terms of the inner product (19), and since 1 a-x(u) = 
1 A(u)e lu ' x , Parseval’s identity gives 


(20) f(x)= E $A{u)\ 2 e iu - x 

u£2nL* 

with convergence for all x. The inversion formula shows that /(it) = |1a(m)| 2 - This 
gives another constraint satisfied by an autocorrelation function, namely that its 
Fourier coefficients are all nonnegative, in other words, / is positive definite. 

We will radialize f by averaging over the sphere or, equivalently, over the or¬ 
thogonal group. In other words, we set 

(21) /(*) = -/ m\\x\\)MO= f f(Tx) dfi(T), 

w (*~ > X ) iS "" 1 /o(R") 

where cu is the surface measure on the unit sphere S'™ -1 C R” and p is the normal¬ 
ized Haar measure on the orthogonal group 0(R n ). 

This function / is radial, i.e., the value of f(x) depends only on ||a:||. (It is not 
periodic any more, nor is it typically the autocorrelation function of a set.) 

Let f l n be the function defined over the nonnegative reals which is such that 

n "(lkll) = 

for all x £ R”. Then, using the inversion formula we get 

f(x) = ^(gn-i) [ _ E f(.u)e zu<M duj(£) = /( M ) fi ™(ll u llll a; ll)- 

^ u(z:2tvL* u(z2ttL* 


f{x) = E«W fi n(t||a;||), 

£>0 


We can rewrite this as 

( 22 ) 
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where n(t) is the sum of f(u) over all u such that ||it|| = t. The sum in (22) is 


absolutely convergent: only countably many of the n(t) coefficients are nonzero, 
they are nonnegative, their sum is 5(A) and the O’s are bounded. 


All constraints in Lemma 3.2 are rotation-invariant (inequality (2) for a fixed 


graph, or inequality (3) for a fixed set is not invariant, but the property that they 
hold for all graphs or sets, respectively, is). As the autocorrelation function of a 
1-avoiding set satisfies these constraints, so does its radialization. We list these 
properties in terms of the function k. Write 5 = 5(A) = /(0). 


(23) 


E t >o K W = <*> 

Et>o = 0, 

Et>o «(*) E* e y(G) fi «(t|M|) < 5a(G) for all graphs G, 

Et>o «(*) E{*, y }e(f) tt n (t\\x - 2 /||) > 5|C| - 1 for all finite C C 


(24) 


We also know that 

«(0) = 5 2 , 

k( t) > 0 for all t > 0. 

In the sequel we will use this (infinite) system of inequalities to get an upper 
bound for 5 and hence for mi(K n ). We will use only a finite number of graphs and 
point configurations, in order to be able to check the properties of the ’’witness” 
function W (f) described below. The above inequalities are necessary but very likely 
not sufficient for a funtion to be the radialization of an autocorrelation function, so 
probably this approach cannot yield the best bound. 

Proposition 3.5. LetS be a finite collection of finite subgraphs of the unit-distance 
graph of R n and let C be a finite collection of finite sets of points in K". Suppose 
that the numbers Vo, v\, wq > 0 for G £ S, and z G > 0 for C £ C are such that 
the function 

(25) W(t) = Uo+UiO„(f)+y^ w G T, Cl n (t\\x\\)-y z c y n n (t\\x~y\\) 

Ges xev(G) cec { x , v }e{ c 2 ) 

satisfies W(0) > 1 and W(t ) > 0 for t > 0. Then toi(R") < 5 where 5 is the 
solution of the equation 

(26) 5 2 = 5 | v 0 + wgoi(G) - y z -c\C\ j 

V Ge5 cec / cec 

Proof. With any function W satisfying W(0) > 1 and W(t) > 0 for t > 0 we have 

(27) 5 2 = k(0) < y<t)w(t). 

t> 0 


If W is in the form (25) then inequalities (|23l) and (27) imply 


(28) 


5 2 < 5 (u 0 + y w G a(G) - y z c \c\ j +yz c . 


Ges 


cec 


cec 


□ 


3.1. Applying Proposition |3.5| for the Euclidean plane. In this section we 
will explain informally how one can look for good collections S and C. Our approach 
for this is experimental. We then give explicit values for v, w, and z in subsection 
3.2 and show how it can be verified that the conditions of Proposition 3.5 are 
satisfied. 
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As we will need to work with the function 122(t), recall that Q n has an expression 
in terms of Bessel functions, namely 

/ 7T\ /2\ (n-2)/2 

^ n(t ) = r^ — J y—J J(n- 2)/2(^) 

for t > 0 and f2 n (0) = 1, where J a is the Bessel function of the first kind with 
parameter a ; this formula was first observed by Schoenberg la¬ 
in trying to apply Proposition 3.5 we can start with S = C = 0. In that case 
W(t) = vq + V\0,2 (t), and the best upper bound that can be achieved this way 
is ~ 0.287, which is slightly worse than the bound of 2/7 « 0.285 coming from the 
Moser spindle (cf. the Introduction). Oliveira and Vallentin [lOj took C = 0 and S 
to be a few equilateral triangles in R 2 (at appropriate positions). This provided an 
upper bound of ss 0.268, which was better than the previously known best upper 
bound of ~ 0.279 due to Szekely [15] . 

A further improvement can be obtained if one takes C = 0 and S to be a few 
congruent copies of the Moser spindle at appropriate positions - as explained here. 
Let G be the Moser spindle as shown in Figure [l] with the lower-left vertex placed 
at the origin. Consider congruent copies of G of the form 

(29) (t,0) + R(9)G, 

where t £ R, 9 e [0, 27 t], and 

f cos 9 — sin 9 


R(9) = 


l sin 9 cos 9 


is a rotation matrix. In order to get a finite number of copies we discretize the 
values of t and 9. Let e = 0.1 and consider —4 < t = je < 4, and 0 < 9 = ke < 2n 
for j, k £Z. As a first step we take S to be all these copies of G. 

Instead of looking for a witness function W ( t ) directly, we will consider the sys¬ 
tem of inequalities (23) and (24). We introduce the normalized variables k(t) = 
n(t)/5, and write the system (23), (24) as a linear programming problem (of infin¬ 
itely many variables k(t)) for maximizing fc(0): 


sup k(0) 


E t > 0 


Et>0 K(t)Q n (t) = 0, 


Et>0 Y2xGV(G) ^n(^ll a 'll) — a (G) 

for Ge5, 

Et>0 «(*) E{x >3/ }e(£) ttn(t\\x - 3/||) > \C\ - 8- 1 

for C G C, 

O 

Al 

-to 

for all t > 0 


(30) 


If for any value of S we obtain sup k(0) = p as the solution of this linear pro¬ 
gramming problem, and p < S, then we should be able to find a witness function 


W(t) in the form (25) by linear duality. In trying to solve this primal problem (30) 
an issue is that we have infinitely many variables k(t). The workaround is to pick 
numbers L > 0 and e > 0 and then only use variables n{t) for t of the form ke < L 
with k > 0 integer. 

By taking L = 200, e = 0.01 and 8 = 0.26305 we get supk(0) »{ = 0.26305. 
Therefore, a witness function W(t) = Vq + rqf l n (t) + w G YY x gV(G) ^n(^ll x ll) 

testifying mi(R") 0.26305 should exist by linear duality (in fact, we must expect 
a little loss because of the discretization in the values of t). Note that this bound 
is already better than that of Oliveira and Vallentin [TO], mi(R n ) < 0.268, using 
equilateral triangles. Also, from the solution of the linear program (30) we can 


see that only a few of the constraints coming from copies of the Moser spindle are 
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actually used; the rest may be discarded. In fact we will only keep three Moser 
spindles, as given in (31) below. 

Now we can try to add inclusion-exclusion constraints. (Note that we have not 
produced any witness function W(t) yet, we have been working with the variables 
k(t). We will keep doing so, and only turn to constructing W(t) at the end.) 

Let k(t) be an optimal solution of the discretized linear program above (con¬ 
taining the constraints coming from the three copies of the Moser spindle given in 
©). Temporarily, we will fix these values k(t) and fix TV > 0. (In our experiments 
the only value that led to improvements was TV = 6 .) We want to find a set C of TV 
points in K 2 which minimizes the left-hand side of the inclusion-exclusion constraint 
in (301. Due to translation invariance, we can assume that the origin belongs to C. 
We can consider the coordinates of the other TV — 1 points as variables, and we 
try to minimize the left-hand side of the constraint using some numerical nonlinear 
optimization method. This will give us a set C of points. To the right hand side 
we substitute the actual value of <5 = 0.26305, and see whether the constraint is 
violated. If it is, then the inclusion of C in C will give us an improved bound on 
mi(R ra ). Namely, we add this new constraint to the linear programming problem 
(30), and find the value of 8 such that the optimum satisfies supre(0) = 5 (note 
here that sup«;( 0 ) depends on 8 because the right hand side of the added con¬ 
straint depends on 8). By linear duality, we can then hope to find a witness W(t) 
which testifies mi(R") 8 (again, a little loss must be expected because of the 
discretization in t). 

This procedure can be iterated: we consider an optimal solution k(i), and find 
another set C of points that minimizes the left-hand side of the inclusion-exclusion 
constraint. And we repeat. After a few iterations, we start to reach the limit of 
this approach. 

Finally we get a dual solution of the resulting linear program (30), which gives 
us the numbers vq, v±, wq and zc to be substituted in Proposition 3.5 These 


numbers will nearly satisfy the conditions of the theorem if our discretization was 
fine enough. After some minor adjustments in these values (incurring a little loss 
in the value of 8) we can actually verify that they satisfy the conditions, and we 
obtain a valid bound for mi(M. n ). 


3.2. Explicit values and verification. We now present candidate values for v, 


w , and 0 for Proposition 3.5 and then show how to modify these values slightly so 
that they satisfy the conditions of the theorem. 

For us, S is composed of three copies of Moser’s spindle, G\, G 2 , and G 3 given 


as in (29), corresponding to the following pairs (f, 6): 


(31) 


(0.4,5.4), (0.6, 5.4), and (0.8, 5.4). 


The collection C we use is given in Table [l] 
We set the values of v, vj, and 0 to: 


( 32 ) 


v 0 2.3022516897351055 
v 1 27.2729338671989154 
Wl 0.2021538298582705 
w 2 0.4311844458316473 
w 3 1.3855315999360112 
Z! 0.2862826361013497 

22 0.7908579212800153 

23 0.9616086568833265 

24 0.2772120180959884 
z 3 0.5311904133936868 
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( 0.781846561681, 0.923983014983) 
(-1.493218191370, 0.715876600816) 
Ci (-0.413794012440, 0.699470697606) 
(-1.195640884520,-0.224511807288) 

( 1.079423910680, -0.016405441239) 

( 0.533352656963, 0.891484779083) 
(-0.611400296245, 0.779442549608) 
C 4 (-0.285536136585, 1.699218505820) 

( 0.251297714466, 0.991412992863) 

( 0.680196169514, 1.919904450610) 

( 0.976422451180, 0.219342709492) 
(-0.896557239530, 0.403173339690) 
C 2 (-0.552919373209, 1.316137405620) 

( 0.051861274857, 0.567345039740) 

( 0.386966521215, 1.028078255990) 

( 0.665906520384, 2.751699047290) 
(-1.358694685180, 1.253666844760) 
C 5 (-0.448270339088, 0.880354589520) 
(-0.943967767036, 0.337966779380) 

( 1.397952081510, 2.088390922180) 

( 0.951509148625, 0.297382071175) 
(-0.856129318724, 0.498561149113) 
C 3 (-0.613074338850, 1.455125134570) 

( 0.035657780606, 0.706699050884) 

( 0.297303535201, 1.051961996200) 



Table 1 . Collection C of point-sets used. Each of the sets also 
contains the origin (0,0), so each set has 6 points. 


Here, w.j is associated with graph Gi in S, and similarly Zi is associated with 
configuration Ci € C. 

Recall that H 2 (t) = Jo(t) and let 

3 5 

= V 0 +ViJ 0 {t) + Y Wi Y MtM)~Y z i Y Mt\\z-y\\)- 

i =1 xeviGi) i= 1 {x,y}e( c 2 -) 

The conditions required of v, w, and z in Proposition |T5] now translate to <p(0) > 1 
and ip(t) > 0 for all t > 0. 

It is easy to check that <p(0) > 1. To show that ip(t) > 0 for all t > 0, the first 
step is to notice that 

lim Jo(t) = 0. 

t—> OO 

This follows from the asymptotic formula for J a for a > 0 (cf. Watson [17], equa¬ 
tion (1) in §7.21). So we see that 

lim ip(t) = v 0 , 

£—>• 00 

and so ip(t) > 0 for all large enough t. 

We need an estimate on how large t has to be chosen and this we can get as 
follows. We have 

dJo(t) _ T 
dt ~ Jl() ' 

Let j 1 < j 2 < J 3 < • • • be the positive zeros of J\. By the above expression for 
the derivative, these are the places of local extrema of Jo- The local extrema of Jo 
decrease in absolute value (cf. Watson [T7], §15.31), that is 

\Jo(jl)\ > I do (j2) | > \Jo(h)\ > ■ ■ 

Hence to find an upper bound on the absolute value of Jo(t) for t > L it is sufficient 
to find the rightmost zero of J\ in the interval [0, L\ and compute Jo at this zero. 
There are procedures to compute the zeros of J\ to any desired precision. 
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Using this idea, we may check that for S and C as we have the absolute value of 

3 5 

ViJ 0 (t) + 

*=1 x £ V ( Gi ) i= 1 { 

for t > 779.8998... (this is the 248th positive zero of J\) is at most vq — 0.05 ss 
2.2522. Consequently ip(t) > 0 for all f > L with L = 780. 

Now we check whether ip(t) > 0 in [0, L\. This will not be the case: since our 
solution has been found numerically via sampling, it will be negative at some points. 
But it will only be slightly negative, and then adding a small number to vq will 
make it nonnegative everywhere. 

Recall that the derivative of Jo is — Ji. Since | Ji(i)| < l/\/2 for all t. > 0 (cf. Wat¬ 
son m, equation (10) in §13.42), we can provide a rough estimate for |<p'(t)|, 
namely 

|<p'(i)| < 75.9547 for allt > 0. 

Then the mean-value theorem implies that 

|<p(*i) - v(te) I < 75.9547|ti - t 2 \ 

for every fi, f 2 > 0. 

If for a prescribed e > 0 we compute ip(t) for all t = ke/76 < L with k > 0 
integer and take the minimum, this gives the minimum of ip in [0, L\ up to an 
additive error of e. 

Taking e = 10 -4 , we obtain the conservative estimate that the minimum of ip 
in [0,L] is at least —0.00011. Adding this to Vo we then get numbers v, w, and z 
that satisfy the conditions of Proposition |3.5[ Solving the quadratic inequality we 
obtain 6 < 0.258795, an upper bound for mi(R 2 ). 

We have attempted to include more constraints from Moser spindles after the 
6 -point configurations were added. This improves the bound slightly, but not much. 
Attempts to add more 6-point configurations have run into numerical trouble, and 
we could not derive a rigorous bound from such trials. They suggest that bet¬ 
ter bounds can be achieved, but we never managed to get below 0.257, which is 
probably the limit of this method. 

The verification procedure we just described was implemented in a Sage m 
script that is available together with the arXiv version of this paper. It is a short 
program that can be easily checked by the reader. Numerical computations are 
still carried out to compute Bessel functions, but due to the simplicity of the code, 
one can have a high degree of confidence on the results obtained. A fully rigorous 
verification procedure would require the use of rational arithmetic and this would 
require Bessel functions to be computed up to good precision using rationals. This 
is not hard to implement, but in our view it is not necessary for the present result. 
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